An approximation to 5' couplings on graphs 
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<^ ■ We discuss a general parametrization for vertices of quantum graphs 

. and show, in particular, how the 5' s and 5' coupling at an n edge 

Qh! vertex can be approximated by means of n + 1 couplings of the 5 type 

■y | provided the latter are properly scaled. 

Quantum graphs became in the last decade a useful and versatile tool to de- 
scribe several classes of physical systems, in particular, various combinations 
.£h ! of quantum wires. There are numerous papers devoted to the subject and 

we restrict ourselves to mentioning the bibliography given in .KS?)!)l IKu04j , 
where also basic concepts of theory are discussed. 

The purpose of this letter is twofold. First of all we want to draw attention 
to useful parametrization of a general coupling at graph vertices whose ad- 
vantages in the present context remained so far unnoticed. Second and more 
important, we address the question of physical meaning of such a coupling 
and suggest an answer illustrating it on a pair of simple nontrivial examples 
of the so-called 5' s and 5' couplings |Ex95| IEx96aj . 

We consider a free spinless particle on a graph, with the Hamiltonian 
which acts as Hipj = —tfi", where ifij denotes the wave function at the jth 
edge. Since early times it has been known that a verte x joining n graph 
edges can be characterized by n 2 real parameters IES89 characterizing the 



1 



boundary condition at the vertex. We use the symbol ^(O) for the column 
vector of the boundary values at the vertex (identified conventionally with 
the origin of the coordinates), and analogously ^'(0) for the vector of the 
derivatives, taken all in the outgoing direction. 

The boundary conditions have to be chosen to make the Hamiltonian 
self-adjoint, or in physical terms, to ensure conservation of the probability 
current at the vertex. A general form of such a coupling was found in [KS99 . 
It is described by a pair of n x n matrices A, B such that rank (A, B) = n 
and AB* is self-adjoint; the boundary values have to satisfy the conditions 

AV(0) + BV'(0) = 0. (1) 

They have an advantage in comparison to earlier parameterizations relating 
^(0) and ^/'(O) by a single matrix, because the latter is typically singular for 
a subset of parameters, albeit a zero-measure one. 

On the other hand, the matrix pair in is non-unique; one would prefer 
to have a condition analogous to -0(0) cos 9 + i/j'(Q)sm9 = is case of a 
single edge end. Such conditions exist, they were obtained independently in 
[FT00, CFTOlj for a generalized point interaction, n = 2, and in jHaOOj for 
any n > 1. It is easy to derive them: the self-adjointness requires vanishing 
of the boundary form, YTj^i^j^'j ~ V^jXO) = 0; which occurs iff the norms 
11^(0) ± i^'(0)||cn with a fixed nonzero £ coincide, so the two vectors must 
be related by an n x n unitary matrix. The length parameter is not important 
because matrices corresponding to two different values are related by 

,_ (£ + £')U + £-£' 
~ {£-£')U + £ + £' 

Thus we set £ = 1, which means a choice of the length scale, and put 

A = U-I, B = i(U + I); (3) 

the edges are obviously fully decoupled at the vertex iff U is diagonal. It is 
easy to check that any such pair satisfies the above quoted requirements from 
KS99J. Conversely, to any A, B with these properties there is a U e U{n) 
and an invertible C such that U = C(A — iB). Indeed, such a U must satisfy 
UU* = C(BB* + AA*)C* since AB* = BA* by assumption. The matrix 
BB* + AA* is strictly positive because its null space is 

kerA* HkerE* = (ran A) 1 - n (ranB) 1 - = (raniU ranS)^ = {0} . (4) 
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In particular, it is Hermitean so C := (BB* + AA*)~ 1 / 2 makes sense, it is 
Hermitean and invertible. 

The parametrization Q simplifies various previous results. For instance, 
the eigenspace of U with eigenvalue —1 gives the projection P = Pi in 
|Ku04| which makes Lemma 4 and the following claims of this paper rather 
transparent. Likewise, the on-shell scattering matrix for a star graph of n 
halflines with the considered coupling equals 

Su{k) ~ (k + i)i + (k-i)u ' (5) 

which makes a discussion of its properties simpler than in Sec. 2 of [KS99 . 

To give an example of the parametrization Q, denote by J the n x n 
matrix whose all entries are equal to one. It is a straightforward exercise to 
check that U = -^j^J — I describes the standard 5 coupling, 

n 

^■(0) = Vfc(0) =: V(0) , j, k = 1, . . . ,n, ^^(0) = «V(0) (6) 

3=1 

with aGl; the case a = corresponds to the "free motion" at the vertex, 
so-called Kirchhoff boundary conditions, while a = oo gives U = —I, the full 
Dirichlet decoupling. In a similar way, U = I — ^z^J describes the singular 
counterpart, so-called 5' s coupling [Ex95, Ex96aj, 

n 

V»5(0) = Vfc(0) =: V^(0) , j,k=l,...,n, ^Vi(0)=^(0) (7) 

3=1 

with Q G R; for /3 = oo we get U — I, the full Neumann decoupling. 

Let us mention another "dual" pair of vertex couplings in which the wave 
functions exhibit permutation symmetry. The more regular one of these is 
the "permuted" 5, or 5 P coupling, given by the boundary conditions 



J2 ^(0) = , ^(0) - ^fc(O) = -(^(0) - ^ fc (0)) , j, k = 1, . . . , n , (8) 
3=i U 

with a G R. It generalizes the 5 S interaction of [TFCO^ and one can check 

easily that the corresponding matrix equals U = Sizi£/ r^J- Its coun- 

terpart is the so-called 5' coupling [Ex95, Ex96aj, 

n 

^(0) = , ^(0) - ^ fc (0) = ^($(0) - ^(0)) , j,k = 1, . . . ,n , (9) 

3=1 U 
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with /3 G E, which corresponds to [/ = —^3^/ + -^z^J- The infinite values 
of the parameters refer again to the Dirichlet and Neumann decoupling of 
the graph edges, respectively. 

Note that in these four examples, the connection condition at the origin 
is totally symmetric with respect to the interchange of edge indices. Conse- 
quently, their U are constructed from symmetric matrices I and J . 

If one wants to continue analysis of such graphs, the first question to be 
answered is about the physical meaning and possible use of the whole family 
of such general couplings. What concerns the second part, a recent inspira- 
tion comes from the domain of quantum computing, where the generalized 
point interactions parameterized by elements of the group U(2) have been 
proposed as an alternative realization of a qubit |CFTQ4j : an extension to 
higher degree vertices opens, of course, interesting possibilities. To make use 
of them, however, one has to understand whether there is a meaningful way 
to "construct" vertices with different couplings. 

The currently available results suggest that this goal cannot be achieved 
in a purely geometrical way, by squeezing a system of branching tubes with 
the same topology as the graph. Several such approximations was analyzed 
recently jKMTl ITTSoTl ISaOTl lEP03| : they all lead either to trivial (Kirchhoff) 
boundary conditions, or to graphs having an extended state Hilbert space 
with extra dimensions due to the vertices. Their common feature was that 
the transverse ground state was a constant function. Hence a nontrivial 
results might be obtained through tubes with Dirichlet boundaries, but this 
problem is open for a long time and notoriously difficult. 

Approximations using potentials scaled in the usual way, i.e. preserving 
their integrals, do yield nontrivial results Ex96b but only for couplings with 
wave functions continuous at the junction, which is just the family (|Sj). This 
is not sufficient and more singular coupling need other means. Our main 
aim here is to explore a natural alternative with approximating interactions 
scaled in a nonlinear way as a generalization of the procedure proposed in 
CS98a, CS98bJ and analyzed from the viewpoint of the convergence topology 
in |AJN00| lEJNZOlj . To keep things simple we will analyze here the 8' s and 5' 
couplings leaving the general case to a subsequent paper. 

Consider first the Hamiltonian Hp on the graph T consisting on n halflines 
coupled at a single vertex by the conditions (JJ|). Consider further the same 
graph with additional vertices of degree two at each arm, all at the same 
distance a > from the common junction. The approximating family will be 
constructed as follows. The operators act, of course, as ipj t— > —ip" at each 
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Figure 1: Scheme of the approximation. For simplicity the graph is featured 
as planar; the vertical bars denote the 5 coupling strength. 

arm; the wave functions satisfy the 5 conditions (JBJ) at the central vertex with 
a coupling parameter a = 6, to be specified later, and another 5 coupling (JBJ), 
this time with the parameter c, at each of the additional vertices - see Fig 1. 
We call such an operator H b,c (a). 

The crucial feature that allows us to simplify the treatment in the present 
situation is a symmetry. Each of the Hamiltonians Hp and H b,c (a) decom- 
poses into a nontrivial part which acts on the one- dimensional subspace of 
7i = ©™ =1 £ 2 (^+) consisting of functions symmetric with respect to permu- 
tations, if)j(x) = ipk{%) for all j, k, and the (n — 1) -dimensional part corre- 
sponding to Dirichlet and Neumann condition at the central vertex for the S 
and 8' s coupling, respectively. Notice that the matrices U corresponding to 
these coupling have each one simple eigenvalue and another one equal to =pl, 
respectively, of multiplicity n — 1. 

To see what the choice of the effective coupling constants 6, c should be, 
let us first modify to our problem the argument of [CS98aj . As we have said, 
in the nontrivial sector all the functions are the same, so we may drop the 
arm index. The boundary values at x = and x = a are related by 

iP(a) = ^(0) + a^'(0) + C(a 2 ) , if/ {a-) = tp'(0+) + 0(a) , (10) 
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V>'(a+) = V'(a-) + c0(a) , V'(0+) = 6^(0) . (11) 

Eliminating "0(0) and ip'(0+) from here, we get in the leading order the 
relation B(a)if)(a) = ^)'(a+), where 

^■■=c+TTab- (12) 
hence the needed limit, /??//(0+) = m/j(0), is achieved as a —>■ 0+ if we choose 

6(a) :=-A, C (a):=-i. (13) 
na z a 

In the orthogonal complement to the permutation-symmetric subspace one 
we can again drop the index, because the operators act in the same way on 
all the linear combinations of X^j=i djifjj(x) which we can choose as the basis 
here, i.e. those satisfying Ylj=i dj — 0. The second one of the conditions (JTTJ) 
is now replaced by ^(0) = 0. Eliminating then the boundary values at x = 
we get in the leading order the relation ip'(a+) — (c+ a~ 1 )V ; («) + 0(a). The 
right-hand side vanishes with the parameter choice (fT3*|) . giving Neumann 
condition, ?//(0+) = 0, in the limit. 

Now we can state and prove our main result. 

Theorem 1 H b ' c (a) Hp as a — >• 0+ in the norm-resolvent sense provided 
the coupling constants b, c are chosen in correspondence with Mfty . 

Proof: By the same symmetry argument as above we can again reduce the 
problem to investigation of a pair of hamine problems. Let us start with the 
one having Dirichlet condition at the origin, so the free Green's function at 
energy k 2 is 

G k (x,y) = S ^^e^>, (14) 

where as usual x< is the smaller one of the values x, y and vice versa. The 
Green's function of the operator with the 5 interaction at x = a is obtained 
easily by Krein's formula AGHH, Appendix A] 

nc( x n ( x . G k (x,a)G k (a,y) 

G k {x,y) = G k {x,y) + — — — — -. (15) 

-c 1 - &k{a,a) 

On the other hand, the Green's function referring to Neumann boundary is 

G?(*,y) = ?^e**>; (16) 
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our aim is to show that the last two converge to each other for some k 2 e C. 
It is convenient to choose k = in with k > 0; we will see below that the 
denominator of the last term at the right-hand side of (JT5j) is then nonzero for 
a small enough. Since the functions involved are uniformly bounded around 
zero, it is sufficient to compute the difference in the case when neither of 
the arguments is smaller than a. For the sake of definiteness suppose that 
a < x < y; then (|T5)l and (jTfij) give 



G c iK (x,y)-Gg(x,y) 



-Ky 



K 



■1 + 



sinh Ka 



-KC~ 



-1 _ q-KX gj^ Ka 



If C 



,-1 



the last term behaves as 1 + 0(a) for a — > 0+, so 
lim G c iK (x,y) = G%(x,y) 

a— >0+ 



(17) 



(18) 



holds for all x, y > 0. 

Consider next the case with the 5 coupling at the origin using the same 
parameter values, namely k = in and a < x < y. We are interested in the 
following two Green's functions, 



Gl(x,y) 
Gi K (x,y) 



-Ky 



' sinh kx + k cosh kx) 



(n sinh + cosh kx) 



(19) 
(20) 



which replace fj!4)l and (|TH|l. respectively, in the present case. The first of 
them determines the full approximating Green's function by Krein's formula, 

G b -(x, y) = G b k (x, y) + ■ (21) 

Using the relations (fTHj) we express the difference 



G b in(x, y) - Gf K (x, y) 



-Ky 



b sinh kx + k cosh kx 



+ 



(b+K) 2 



' sinh kx + k cosh kx) 



6+ K 

n sinh kx + (3k cosh kx 



Ka 



— (b sinh kx + k cosh kx) 



(22) 



The first term in the bracket tends to sinh hZOC clS (X 
is independent of a, so their sum in the limit gives 

(3Ke~ KX 
n + (3k 



n + (3k 
0+, while the third one 

(23) 
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Next we take the middle term without the factor e~ KX and expand the nu- 
merator and denominator to the second power in a; this gives its limit which 
differs just by the sign from (J2"3|) . in other words 

\im Gi c (x,y) = Gi(x,y) (24) 

a— >0+ 

holds again for all x, y > 0. To conclude the proof we have just to realize 
that as functions of x, y the differences (JT7J) and (|2*2*jl decay exponentially, so 
the corresponding resolvent differences converge to zero even in the Hilbert- 
Schmidt norm. □ 

Let us add that the proven result opens way to approximation of Hp by 
Hamiltonians with more regular interactions. We have mentioned that the 
central S in H b,c (a) can be approximated by a family of potentials scaled in 
the usual way, the same is true for the 5 interactions at the points Xj = a. As 
in the related problem discussed in |EJNZ01j . the question then is how fast 
have these approximating potentials to shrink with respect to a. 

Consider finally the case of 5' coupling, i.e. the Hamiltonian Hp on our 
star graph with the boundary conditions Q at the vertex. The approximat- 
ing family will be constructed in a similar way as above; the difference is that 
now the wave functions will satisfy the 5 P conditions (jHJ) at the central vertex 
with a coupling parameter a = 6, to be specified. The rest is the same, there 
is another 5 coupling © with the parameter called again c at each of the 
additional vertices; we denote such an operator H b,c (a). 

To realize that this problem can be again reduced to a one-dimensional 
analysis, denote e := e 2m / n and introduce 

g r := {(iPix),^^), . . . ^^-^(x)) : ^GL 2 (R + )} . (25) 

The graph state Hilbert space can be written as H = ®"Zq Q r - Indeed, any 
vector of Ti is a unique linear combination of the vectors from Q r , because 
the determinant of the corresponding linear system is the Vandermond de- 
terminant of 1, e, ... , e n_1 , which is nonzero because the latter are mutually 
different. It is straightforward to see that the subspaces Q r are invariant un- 
der the Hamiltonians in question; the 5 p and 5' couplings acts "trivially" at 
the origin corresponding to Dirichlet and Neumann condition, respectively, 
while on each of the subspaces Q r , r = 1, . . . , n, these boundary conditions 
are replaced by ^'(0) = ~ip(0) and ^(0) = ^'(0). Thus we can choose 

6(o):=-4, c(o):=--, (26) 
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and repeat the above considerations, arriving at the following conclusion. 

Theorem 2 H b,c (a) — ► Hp holds in the norm-resolvent sense as a —>■ 0+ if 
the coupling constant families b, c are given by Wty . 
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